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ABSTRACT 

The formation of CH + in the interstellar medium has long been an outstanding problem in chemical 
models. In order to probe the physical conditions of the ISM in which CH + forms, we propose the use 
of CHj observations. The pathway to forming CH^ begins with CH + , and a steady state analysis 
of CH^ and the reaction intermediary CH^ results in a relationship between the CH + and CH^ 
abundances. This relationship depends on the molecular hydrogen fraction, /h 2 , and gas temperature, 
T, so observations of CH + and CH^ can be used to infer the properties of the gas in which both species 
reside. We present observations of both molecules along the diffuse cloud sight line toward Cyg OB2 
No. 12. Using our computed column densities and upper limits, we put constraints on the /h 2 vs. 
T parameter space in which CH + and CH^ form. We find that average, static, diffuse molecular 
cloud conditions (i.e. /h 2 > 0.2, T ~ 60 K) are excluded by our analysis. However, current theory 
suggests that non-equilibrium effects drive the reaction C + + H2 — >• CH + + H, endothermic by 4640 
K. If we consider a higher effective temperature due to collisions between neutrals and accelerated 
ions, the CHjj~ partition function predicts that the overall population will be spread out into several 
excited rotational levels. As a result, observations of more CH^" transitions with higher signal-to-noise 
ratios are necessary to place any constraints on models where magnetic acceleration of ions drives the 
formation of CH + . 
Subject headings: astrochemistry 



1. INTRODUCTION 



1.1. 



Background 

Although CH + wa s first discovered in interstellar space 
nearl y 70 years ago (|Dunhaidll937tlDouglas fc Herzberd 
119411 ). the mechanism by which this simple molecule 
forms has remained elusive. This is because many the- 
oretical chemical models have been unable to reproduce 
the large observed abundance of CH + in diffuse cloud 
sight lines. At present, it is thought that CH + is primar- 
ily formed by the reaction 



Crawford 1995), turbulent dissipation and Alfven waves 
in diffuse clouds remain viable mechanisms by which the 
rate of reaction ([lj may be increased, and there is no 
clear reason to favor one theory over the other at present. 



1.2. 



CH^~ Chemistry 



C H 



H 2 



CH 



H. 



(1) 



However, this reaction is highly endothermic (Ap~ [ = 1.0 x 
10" 10 exp(-4640/T) cm 3 s -1 ; iFederman eFaLl 119961 ). 
such that non-equilibrium chemistry must be invoked in 
order for it to proceed rapidly enough to produce the 
observed amounts of CH + , without vastly overproducing 
observed abundances of OH. 

Over the past several decades, various theories have 
been proposed to account for an increased rate of reaction 
(Ht, including, but not l imited to, neutral shocks (e.g. 
lElitzur fc Watsonl I1978L I1980D. magneto hvdrodvnamic 



To place constraints on the physical conditions of the 
interstellar medium (ISM) where CH + forms — and po- 
tentially discriminate between the various proposed for- 
mation mechanisms — we have investigated the reaction 
network linking CH + and CH^", and made/obtained ob- 
servations searching for both species. In diffuse molec- 
ular clouds the carbon chemistry is either initiated by 
reaction ([TJ , or by the radiative association reaction 



C^ 



H, 



CH+ + hv, 



(2) 



which "bypasses" CH + production. The branching frac- 
tion between these reactions is temperature dependent, 
and reaction ([1]) will dominate when collision tempera- 
tures are greater than about 400 K. Assuming CH + is 
formed, it will react with molecular hydrogen, 



(MHD) shocks (e.g. iDraine fc Katzl ll986[). Alfven waves 
(IFederman et aLl 11996ft. and turbulent dissipation (e.g. 
Dulev et al.l Il992t iGodard et al.l 120091: iPan fc Padoanl 
2009). While both neutral and MHD shocks seem to 



CH 



H, 



cm 



and eventually form CH3 via 



have been ruled out by observations (jGredel et alJll993l 
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In addition to reaction ((4]) though, the CH^~ intermedi- 
ary is also destroyed by dissociative recombination with 
electrons 



cm 



e — > products, 



(5) 



thus decreasing the CHg production rate. Dissociative 
recombination with electrons also happens to be the pri- 



2 



mary process by which CHJj is destroyed: 

CH^ + e" -t products (6) 

(there are multiple product channels for reactions (|5j) 
and ([6|), but our analysis only depends on the over- 
all dissociative recombination rates). While CH^" will 
be formed starting from both reactions ((2]) and ((3]), 
the latter process will dominate in clouds contain- 
ing CH+ (/^(CH+)/A fop;(C + ) > 100 whe re = 
1.2 x 10~ 9 cm 3 s" 1 pdcEwan et"aLl fl999t) ^ = 
4 x lQ- 16 (T/300) - - 2 cm 3 s" 1 ([Herbstl IIMI fl985h . 
x(CH+) ~ 10~ 8 (ISheffer et all I2008I) . and x(C+) ~ 
1.4 x 10" 4 (jCardelli et al.lll996h ). Because we are study- 
ing cloud components with CH + , we omit reaction ([2]) 
from our analysis. Note, though, that if CH^" were ob- 
served in a cloud component lacking CH + , the following 
analysis would not be applicable. 

Assuming steady state for both CH^~ and CH^ , we can 
derive a relation between the concentrations of CHjj" and 
CH + , given by 

"(CH+) = fafa 

n(CH+) 2^p e (/ H2 + 2x e kg}^' [) 

Here, the k^s are the rate coefficients for reaction i, 
/h 2 = 2n(H2)/nH is the molecular hydrogen fraction 
(where hh = n(H) + 2ri(H2)), and x e is the electron 
fraction (x e = n e /nn). The rate coefficients used in this 
study are Agj = 1.6 x 10" 9 cm 3 s" 1 ([Smith fc Adamsl 
19771 ), /j5|= 6.4xl0" 7 (r/300)" - 6 cm 3 s" 1 ([Larson et afl 
1991, and fm\= 3. 5 x 10" 7 (r/300)- - 5 cm 3 s " 1 (IMitchelll 
19901) . However, iSheehan fc St. -Maurice! ([2004J) note 
that to] was determined using vibrationally excited CH^ . 
The rate coefficient for dissociative recombination from 
the ground vibrational state — the only state likely to 
be populated in diffuse cloud conditions — may differ 
from the above experimental value by about a factor of 
3. Because of this possibility, during our analysis we also 
consider using 3^jg]and ^gj/3 in equation ([7]). 

Aside from the electron fraction — estimated to be 
x e ~ 1.4 x 10~ 4 in diffuse clouds via C + observations 
(|Cardelli et al.l [19961) — equation ([7|) is dependent only 
on the molecular hydrogen fraction and the temperature 
(through the rate coefficients ^ and t^. As a result, 
the ratio between the abundances of CHjj" and CH + is 
governed by the interstellar physical parameters /h 2 and 
T. Observations of CH + and CH^ can thus be used 
to infer the physical conditions in the ISM where these 
molecules reside. 

2. OBSERVATIONS & DATA REDUCTION 

2.1. The Cyg OB2 No. 12 Sight Line 

We selected the Cyg OB2 No. 12 sight line for this 
study becaus e it combines a large extinction (Ay ~ 10; 
ISchultd Il958f ) and thus gas column, with a background 
star that is bright enough at both optical and infrared 
wavelengths such that the observations described be- 
low are possible. However, there is some contention 
as to whether the sight line mainly probes diffuse or 
dense material, environments which have different elec- 
tron fractions. Rotational excitation analyses have been 



performed using C2 observations, and the resulting in- 
ferred hydrogen number densities are between 200 cm" 3 
and 4 00 cm" 3 (|Gredel et al.ll200H iSonnentrucker et al.l 
120071) . typical of diffuse cl ouds. Also, the non-d etection 
of H 2 and C0 2 ices by IWhittet et ail ([1997D led to 
the conclusion that there was a lack of dense molecu- 
lar gas along the sight line. On the other hand , ra- 
dio ob servations of HC O + ([Scappini et al.l 120001) and 
13 CO (jCasu et al.l 120051 ) suggest that portions of the 
sight line may pass through denser c lumps as pro- 
posed b y iCecchi-Pestellini fc DalgarncJ (120001) . Addi- 
tionally, jCecch^esteffin^^^algan^ that 
dense clumps can still be present even with the afore- 
mentioned C2 analyses. Because we are examining the 
chemistry related to CH + , a species thought to reside 
mainly in diffuse gas with 10 cm" 3 < nn < 300 cm" 3 
(|Pan et al.ll2005D . we adopt the diffuse molecular cloud 
model for our analysis. However, in Section 4 we also 
consider a reduced electron fraction in order to investi- 
gate the effects of varying this parameter. 

2.2. CH+ 

Observations toward Cyg OB2 No. 12 and the at- 
mospheric standard a Cyg we re made on 1999 Nov 19 
using the CGS4 spectrometer ([Mountain et al.l [19901 ) at 
UKIRT. The spectrometer was used with its echelle grat- 
ing, 0.6" wide slit, and long camera to yield a resolv- 
ing power of about 37,000, and a circular variable filter 
(CVF) was employed to select the correct order. Spec- 
tra were centered to cover the r i?(l,0), r R(l,l), and 
r R(2, 2) ro- vibrational transitions of CH^ at 3.18523 /im, 
3.19578 /x m, and 3.18777 /an (vacuum wavelengths), re- 
spectively (Crofton et al.11985?) . Stars were nodded along 
the slit in an ABBA pattern with total integration times 
of 864 s for Cyg OB2 No. 12 and 672 s for a Cyg. 

Individual frames were processed using Starlink's 
ORAC-DR pipelin^E specifically designed to handle 
UKIRT data. Spectra were then extracted from the re- 
sultant group frame (the combination of all integrations 
for a given target) using NOAO's IRAF packag^l, and 
imported to IGOR Pro3 where we hav e macros set up 
to complete the reduction (|McCallll2001[ ). Here, artifacts 
from the data acquisition methods were removed, spectra 
were wavelength calibrated using telluric lines, and the 
science target spectra were divided by the telluric stan- 
dard spectra to remove atmospheric absorption features. 
The resulting normalized spectrum of Cyg OB2 No. 12 
is shown in Figure [TJ 

2.3. CH+ and CH 
The A-A(O-O) and A-X(l-0) transitions of 
CH+ (4232.548 A and 3957.692 A, respectively; 
ICarrington fc Ramsayl Il982l ) and the A-X(0-0) transi- 
tion of CH (4300.308 A; lZachwieialll995D . were observed 
simultaneously u sing the HIRES instrument at Keck 
([Vogt et all I1994D by G. Blake and collaborators. Ob- 
servations toward Cyg OB2 No. 12 were performed on 
1999 Jun 17 and 18 with a slit width of 0.57", producing 

5 See web site at ht tp:/ /www.oracdr.org/[ 

6 See web site at http://iraf.noao.edu/. 

7 See web site at http://www.wavemetrics.com/. 
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3.190 3.192 3.194 3.196 

Wavelength (fim) 

Fig. 1.— This spectrum of Cyg OB2 No. 12 covers the r R(l, 0), 
r R(l, 1), and r R{2, 2) transitions of CH+. The three vertical lines 
underneath each transition label show the expected positions for 
absorption due to various cloud components along this line of sight. 
The > 2% fluctuations near 3.185 jim, 3.1965 /im, and between 
3.191 (im and 3.194 fim are due to strong telluric absorption fea- 
tures which could not be completely removed through standard star 
division. The feature near 3.189 fim is an instrumental artifact. 

a resolving power of about 70,000. The total integration 
time on source was 8400 s. 

Using IRAF, the profile of the overscan region was sub- 
tracted from each image, and the frames were averaged 
together with the cosmic ray reject option enabled. The 
relevant orders (10, 16, & 17) were extracted using apall 
with background subtraction enabled^, and normalized 
using the continuum routine. These astronomical spectra 
were then imported to IGOR PRO, along with Thorium- 
Argon arc lamp spectra used for wavelength calibration. 
After calibration, observed wavelengths were converted 
to LSR velocities, and the resulting spectra are shown in 
Figured 

3. RESULTS 



3.1. 



CH+ 



The spectrum in Figure [T] shows no indication of ab- 
sorption from any of the CH^ transitions considered. 
Upper limits on the equivalent widths were computed 
using W\ < crApi x ^/A/pi x , where a is the standard de- 
viation on the continuum of the spectrum, A prx is the 
wavelength per pixel, and A/" p i x is the number of pixels 
expected in an absorption feature given a full width at 
half-maximum (FWHM) of 12.4 km s _1 (average FWHM 
measured for the CH + and CH lines). The upper limits 
on column densities were then calculated using the stan- 
dard relation between equivalent width and column. The 
results of this analysis are shown in Tabic Q] 

3.2. CH+ and CH 

The CH + and CH spectra in Figure [2] show broad 
and deep absorption profiles. While it is known that 
there are multiple velocity components toward Cyg OB2 

8 The initial background subtraction for order 10, which con- 
tains the 3957.692 A line of CH + , resulted in intensities below 
zero at the centers of Ca II lines near 3934 and 3968 A. This is 
most likely due to the low flux and low signal-to-noise ratio of this 
order. A curve of growth analysis (see Section 3.2) using the equiv- 
alent width of the 3957.692 A line from this uncorrected spectrum 
resulted in no value of the Doppler parameter for which the column 
densities determined from the A— X(0-0) and A— X(l-0) absorption 
lines agreed. To remedy these unphysical results, we shifted the 
spectrum of order 10 so that the deeper Ca II line had zero intensity. 




LSR Velocity (km s" ) 

Fig. 2.— Spectra of Cyg OB2 No. 12 showing the A-X(0-0) 
transition of CH (top spectrum) and the A-X(0-0) and A— X(l- 
0) transitions of CH + (middle and bottom spectra, respectively). 
The vertical dotted lines mark cloud compon ents at 6.9 km s — 1 
and 12.3 km s" 1 seen by IMcCall eTaTI 120021 1 in CN, C 2 , 12 CO, 
Hg~, and K I. The dashed line marks a 1 km s -1 component only 
seen in Hi" and K I. 



No. 12, the low signal-to-noise ratios of the CH and 
CH + spectra make fitting multiple components to the 
observed absorption features highly uncertain. Because 
there is no unique way to fit any of the absorption fea- 
tures with 3 separate components, we simply found the 
equivalent widths of each feature as a whole. Using 
these equivalent widths, we first computed lower limits 
to the column densities of both species by assuming that 
the lines are optically thin. As this assumption is most 
likely incorrect, we also computed column densities us- 
ing a curve of growth analysis. Both CH + lines arise 
from the same lower energy state, and so should result 
in the same column density. Plotting the column den- 
sity vs. the Doppler parameter, 6, for both lines, we 
found that the curves intersect at b = 3.5^4 km s _1 and 
N(CH+) = 2.2lg j x 10 14 cm" 2 within the uncertainties 
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of the equivalent widths (the positive uncertainty asso- 
ciated with b was computed from the average FWHM of 
the CH + lines and the relation b = Awfwhm/ (2-v/ln(2))). 
Assuming that b = 3.5^ km s _1 is also applicable to CH, 
we computed 7V(CH) = 2.8^ x 10 14 cm" 2 . All of these 
results, as well as the various line parameters, are shown 
in Tabled 

3.3. Interstellar Conditions 

Using both CH + and CH^ observations, and as- 
suming that the two species have the same distribu- 
tion along the line of sight (i.e. N (CH+ ) / TV (CH+ ) = 
n(CH^)/n(CH + )), equation (JTJ) can be used to constrain 
the temperature and molecular fraction of the gas where 
these species reside. However, equation (J7J requires 
the total abundance of CH;j~, whereas we only have up- 
per limits for the column densities of three particular 
states. As a result, the equation must be recast such 
that the left hand side is the ratio between the column 
density of an individual state of CHjj~ with the total col- 
umn of CH + . This is achieved by using the partition 
funct ion for CHj~ (compu ted using molecular co nstants 
from ICrofton et all 119881 and IJagod et all 1 19941) to cal- 
culate the fractional population, Pj,k(T), in each indi- 
vidual state as a function of temperature. The defini- 
tion AT(CHt )Pj,k(T) = N(J, K) is then substituted into 
equation giving 

iV(CH+) JKV >21 We {fK 2 +2x e l^l^§- y> 

Using the observed CH + column and the individual 3c 
upper limits for CH;j~, we can put three separate up- 
per limits on the left hand side of the equation. The 
(/h 2 ,T) parameter space can then be explored to de- 
termine which molecular fraction-temperature combina- 
tions are excluded by our observations. This analysis is 
shown in Figure |3] 

4. DISCUSSION 

It is quite clear from Figure [3] that our analysis us- 
ing CH + and CHj observations excludes a large portion 
of the (/h 2 , T) parameter space. Indeed, typical diffuse 
molecular cloud conditions (fu 2 ^ 0.2, T ~ 60 K) are ex- 
cluded by our upper limits to the column densities of all 
three observed CH^ transitions. While this may not be 
especially surprising — equilibrium chemistry in diffuse 
clouds is unable to account for the large abundances of 
CH + — our observations add a new, independent method 
for determining that CH + cannot form under such con- 
ditions. 

However, as mentioned in Section 1.2, the dissociative 
recombination rate coefficient for CH^~ with electrons, 
fej, is most likely uncertain by a factor of about 3. If 
we allow variations between 3 fen and fej]/3 during our 
analysis using JV(1,0) with equation ((8|), the dashed and 
dotted contours in Figure [3] result, respectively. With 
fej/3 essentially every temperature is excluded when the 
molecular fraction is greater than 0.1. For 3fep. high 
molecular fractions are allowed as the temperature in- 
creases beyond a few hundred Kelvin. Even with this 
higher rate coefficient though, most typical diffuse molec- 
ular clouds conditions are still excluded by this analysis. 
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Fig. 3. — Plot of the (/h 2 >^) parameter space. The 3 solid 
curves represent contours corresponding to the values of /jj 2 and 
T which reproduce the ratios N(l, 0)/N(CH+), N(l, 1)/JV(CH+), 
and N(2, 2)/W(CH+)_( where the N(J, K) values are 3tr upper lim- 
its) using equation J5J, and are labeled accordingly. Shaded re- 
gions are excluded by our analysis, and progressively darker shad- 
ing indicates regions excluded by more than one observation. The 
dashed and dotted curves show the parameter space excluded by 
the ./V(1,0) upper limit using 3/jg]and ^g]/3, respectively. 

While we have assumed diffuse cloud conditions 
throughout our analysis and given reasons for doing so, 
we feel it is prudent to revisit the dense clump scenario 
described in Section 2.1. For the sake of completeness, 
we examine the effects that a lower electron fraction 
(due to denser material) would have on our analysis. If 
we decrease the electron fraction by a factor of 10 to 
x e = 1.4 x 10~ 5 (still m uch higher than x e ~ 4 x 10~ 8 
assumed in dense clouds (jWoodall et al.l l2007)). then in- 
terstellar clouds with /h 2 > 0.04 and 10 K < T < 1000 K 
are excluded. This demonstrates that as x e decreases, re- 
gions with even smaller molecular hydrogen fractions are 
ruled out by our analysis. 

A potentially larger uncertainty in our analysis is that 
we only consider a single kinetic temperature of the gas, 
and do not account for any non-equilibrium effects. If 
ions are magnetically accelerated relative to neutrals, our 
analysis changes drastically. Assuming an effective tem- 
perature for collisions can b e described as (|Flower et al.l 
[19851: iFedeFman et al.lU996l ) 

Tea = T khl + ^-(Av) 2 , (9) 
Sk B 

where /x is the reduced mass of collision partners, ks 
is Boltzmann's constant, and Av is the turbulent veloc - 
ity (about 3.3 km s _1 accoring to She ffer et al.l l2008h . 
the effective temperatures for collisions of CHj^ with H2 
and H are about 800 K and 400 K higher, respectively, 
than the kinetic temperature of the gas. These higher 
temperatures are important because it is generally as- 
sumed that the fractional populations of CH^" states are 
"thermalized" by collisions with H and H2. If collision 
temperatures are hundreds of Kelvin instead of ~60 K, 
then the expected populations of the (1,0), (1,1), and 
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(2,2) states we observed are significantly decreased. In 
this case, our analysis is unable to exclude any significant 
portion of the (/h 2 ,T) parameter space. 

To use our method of analysis assuming high effective 
temperatures due to accelerated ions, new observations 
must be made. Such observations should cover the 3 
transitions we have already examined, search for absorp- 
tion due to higher energy states of CK^, and obtain a 
higher signal-to-noise ratio in all cases. If absorption 
from CHg can be detected, instead of only excluding 
portions of (/h 2 ,T) space, we will be able to constrain 
the conditions where CH + does form. 

5. CONCLUSIONS 

We have used observations of CH^ and CH+ toward 
the diffuse molecular cloud sight line Cyg OB2 No. 12 
in an attempt to constrain the physical conditions of the 
ISM where these species reside. If we assume equilibrium 
chemistry, our analysis excludes the portion of (/h 2 ,T) 
parameter space corresponding to average diffuse molec- 
ular clouds. This acts as an important independent check 
on previous studies which also concluded that the slow 
rate of CH + formation ruled out equilibrium chemistry 
in such environments. However, our analysis cannot ex- 
clude non-equilibrium effects in diffuse clouds, and nei- 
ther supports nor refutes any of the current proposed the- 
ories for CH+ formation. To do so, new observations of 
CHjj" with signal-to-noise ratios much higher than those 
obtained in this study (~ 600) would be required. Also, 



it would be highly advantageous if the dissociative re- 
combination rate coefficient for CHj were measured at 
low vibrational temperatures, i.e. for the ground vibra- 
tional state. 
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TABLE 1 
CHg" Absorption Line Parameters 





M 2 


W x 




Transition 


(D 2 ) 


(mA) 


(10 12 cm" 2 ) 


r K(l,0) 


.0313 


< 1.34 


< 3.23 


r _R(l, 1) 


.0313 


< 1.21 


< 2.91 


r R(2,2) 


.0313 


< 1.10 


< 2.65 



Note. — Values for W\ and JVi cvc i are lcr 
upper limits. Transition dipole moments were 
calculated using Honl-London factors and 
the tr ansition moment given bv lPracna et all 
(1993) 
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TABLE 2 

CH and CH+ Absorption Line Parameters 



Species 


Transition 


/ 


W x 
(mA) 


FWHM 

(km s— 1 ) 


^thin 

(10 14 cm- 2 ) 


N 

(10 14 cm" 2 ) 


b 

(km s -1 ) 


CH+ 
CH+ 
CH 


A-X(O-O) 
A-X(l-O) 
A-X(O-O) 


0.00545 a 
0.00331 a 
0.00506 b 


103 ±6 
70 ±8 
113 ±4 


13.7 
11.3 
12.2 


> 1.2 

> 1.5 

> 1.4 


9 9+1.5 

*-*-o.e 

9 9+1.5 
^-0.6 

2.8l 4 


3.5 +4 
3.5 +4 
3.5 +4 



Note. — Lower limits (JVthin) were computed under the assumption that the gas is optically 
thin. As this is most likely not the case, we performed a curve of growth analysis for both CH + lines 
(which arise from the same lower state) in order to find the b- value where both lines predict the same 
CH+ column density. This occurs at b = 3.5^ 4 km s~\ where JV(CH+) = 2.21q J X 10 14 cm" 2 . 



a Oscillator strength from Larsson & Sicgbahn (1983a). 
b Oscillator strength from Larsson & Sicgbahn (1983b). 



